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Abstract 

We review the notion of differential Fay identities and demonstrate, 
through case studies, its new role in integrable hierarchies of the KP 
type. These identities are known to be a convenient tool for deriving 
dispersionless Hirota equations. We show that differential (or, in the case 
of the Toda hierarchy, difference) Fay identities play a more fundamental 
role. Namely, they are nothing but a generating functional expression of 
the full set of auxiliary linear equations, hence substantially equivalent 
to the integrable hierarchies themselves. These results are illustrated for 
the KP, Toda, BKP and DKP hierarchies. As a byproduct, we point out 
some new features of the DKP hierarchy and its dispersionless limit. 
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1 Introduction 



The simplest, but most important example of dispersionless limit can be seen 
in the KdV equation 

dtu + GudxU + e'^d^u = 

with a small dispersion parameter in front of the dispersion term. By setting 
e = 0, the third order PDE reduces to the first order PDE (the dispersionless 
KdV equation) 

dtU + GudxU = 0. 

This naive procedure changes the properties of solutions drastically. For exam- 
ple, since soliton solutions can exist in the balance of nonlinearity and dispersion, 
the dispersionless KdV equation no longer has soliton-like solutions. Nonlinear- 
ity without dispersion leads to the so called gradient catastrophe, namely, a 
regular initial value develops singularity ("shock") with \dxu\ = oo. This means 
that the appoximation by the dispersionless equation breaks down at that stage. 
In the presence of a small but nonzero dispersion term, this singularity is reg- 
ularized to form an oscillatory region called "dispersive shock" (see the review 
by Lax, Levermore and Venakides [T]). Nevertheless the dispersionless KdV 
equation is a good approximation until the gradient catastrophe takes place. 
Moreover, this equation inherits many aspects of integrability of the KdV equa- 
tion such as the existence of a Lax formalism, an infinite number of commuting 
flows, etc. We can consider a dispersionless limit of the KP equation 

dyU = dx{dtu + udxU + e^d^u) 

in a similar manner. 

Another significant example is the two-dimensional Toda equation. This is 
originally a system of PDEs 

dxdy(l>s + e'^=+i-'^= - e-^^-^^-i = 

for an infinite number of fields (jfs — (j)six,y) indexed by an integer s S Z. x 
and y are light cone coordinates of the two-dimensional Minkowski spacetime. 
If s is a coordinate on a one-dimensional lattice, and (f>s is interpreted as a field 
(/)(s) = (j3{s,x,y) on a partially discretized three-dimensional spacetime, we can 
consider a continuum limit as the lattice spacing e tends to 0. To this end, we 
restart from the rescaled Toda equation 

edxdy<l>{s) + e("*(^+')-^(^))/^ - eW(^)-^(^-^»/^ = 

and let e ^ 0. This yields the dispersionless Toda equation 

dxdy^ + dsie^^'f') = 0. 
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Remarkably, this equation is also known in general relativity as the Boyer-Finley 
equation, which was discovered in a classification of selfdual spacetimes with a 
Killing symmetry This equation, too, is a kind of integrable system. 

These rather naive constructions of dispersionless limit can be reformulated 
in a more systematic way in the framework of the KdV, KP and Toda hierarchies. 
The outcome are the the dispersionless KdV, KP and Toda hierarchies. These 
dispersionless integrable hierarchies share many properties with the original 
dispersive systems. We refer details to the review by Takasaki and Takebe |4] 
and recall a few essential features. 

Firstly, the dispersionless limit turns out to be an analogue of "quasi-classical 
limit" in quantum mechanics. The small parameter e amounts to the Planck 
constant h. This analogy (quantum-classical correspondence) becomes particu- 
larly fruitful when one considers the Lax formalism in the dispersionless limit. In 
the dispersionless KP and Toda hierarchies, commutators of (pseudo)differential 
operators in the Lax and Zakharov-Shabat equations are replaced by Poisson 
brackets, and the role of auxiliary linear equations are played by Hamilton- 
Jacobi equations. Actually, these ideas are rather old and can be found in 
studies of the Benney hierarchy in the early eighties O [H [7] . 

Secondly, one can formulate the quasi-classical limit in terms of the tau 
function as well. This is based on ideas borrowed from random matrices, topo- 
logical field theory and string theory [51 [SI [10] • Although the Hirota equations 
themselves do not survive the dispersionless limit, one can find an alternative 
framework, namely, "dispersionless Hirota equations" [TTl[4l[T2]. Actually, we 
can see a prototype of dispersionless Hirota equations in the aforementioned 
dispersionless KdV and Toda equations. In the case of the dispersionless KdV 
equation, we can convert the equation to 

dtd^F + 3idlFf =0 

by changing variables as m = 2d^F and integrating the equation once with 
respect to x. In the case of the dispersionless Toda equation, we can similarly 
derive the equation 

d.^dyF + e^'^ = 

by subtstituting (f) — dgF. The new dependent variable F may be thought of 
as a counterpart of the tau function in the dispersionless systems. Note that, 
unlike Hirota equations in the dispersive case, the equations for the F function 
are no longer bilinear. 

Dispersionless Hirota equations have played a central role in recent studies 
on dispersionless integrable hierarchies. They were applied to various problems 
of mathematical physics such as interface dynamics ,13] , associativity equations 
[13] and string field theory [TSl [HI [13 [H] • Alongside these apphcations, math- 
ematical aspects of dispersionless Hirota equations were also studied in detail 
[ll[20l[2Tl[22l[2l[24l[2Sl[26|. It has been estabhshed through these studies that 
dispersionless Hirota equations can be a fundamental language for dispersionless 
integrable hierarchies. 
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By the way, the notion of dispersionless Hirota equations was first discovered 
(in the case of the KP hierarchy [H Appendix B]) as dispersionless hmit of the 
so called "differential Fay identity" . In other words, this identity (introduced 
by Adler and van Moerbeke in a quite different context) is a dispersive 
counterpart of dispersionless Hirota equations. This raises a natural question: 
Can the notion of differential Fay identity play the same fundamental role as 
that of dispersionless Hirota equations? 

This paper presents several case studies on this question. In the case of 
the KP hierarchy, this question was already answered affirmatively when the 
notion of dispersionless Hirota equations was first proposed [1] . It was proven 
therein that the differential Fay identity, which is a consequence of the KP 
hierarchy, is actually equivalent to the KP hierarchy itself. We shall review this 
result, and point out that the most essential part of this result is the fact that 
the differential Fay identity is a generating functional expression of an infinite 
number of auxiliary linear equations. We shall show that the same point of view 
is valid for some other cases as well. 

This paper is organized as follows. Sections 2 and 3 are devoted to the most 
fundamental cases, namely, the KP and Toda hierarchies. All results presented 
here are well known and can be found in the literature. The case of the KP 
hierarchy is reviewed in detail (Section 2) as a prototype of the subsequent 
cases. In the case of the Toda hierarchy (Section 3), difference analogues of 
differential Fay identities show up. Sections 4 and 5 deal with the BKP and 
DKP hierarchies, which are relatives of the KP hierarchy. Most results on the 
BKP hierarchy (Section 4) are already published except for the relation between 
the differential Fay identities and the auxiliary linear problem. The results on 
the DKP hierarchy (Section 5) are mostly new. In Section 6 we show our 
conclusion along with several remarks. 

2 KP hierarchy and differential Fay identity 
2.1 Bilinear equations for tau function 

Let t — {ti,t2, ■ ■ .) denote the set of time variables of the KP hierarchy. The 
first one ti is identified with the spatial variable x in the Lax formalism. 

The tau function t = T{t) of the KP hierarchy [28l [29l |30l |31] satisfies 
an infinite number of bilinear equations. These equations can be written in a 
compact form as 

{Df ~ AD1D3 + 3DI)t ■ T = 0, .... (2.1) 
with Hirota's notation 

P{Di,D2, ...)t-t= P{d[ - ai, - 92, . . .)^(ti,4, • ■ ■)r{tut2, . . .)\,,^, , 

where i9„ and d'^ denote the derivatives (9„ = d/dtn, d',^ — djdt'^ in the compo- 
nent of t and t' = (t']^, • • ■)• 
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Date, Jimbo, Kashiwara and Miwa discovered that these Hirota equations 
can be encoded to a single equation of the form 

^e«(*'-*-)r(t'-ri]Mt + [z-i])^0, (2.2) 

where £,{t,z) and [a] are the standard notions 

^it,z) = ^tnz"^, [a] = ^'a, Y> ■ • ■ ' • ■ •) . 

n=l ^ / 

and f is the contour integral along a sufhciently large circle |z| = i? (or a formal 
algebraic operator extracting the coefficient of z^^ of a Laurent series). (|2.2p 
is understood to hold for arbitrary values of t' and t. More precisely, (|2.2p is a 
generating functional expression of an infinite number of equations that can be 
obtained by Taylor expansion at t' = t. 

One can convert (|2.2p into a form that is more directly related to Hirota 
equations. This is achieved by introducing a new set of variables a = (ai, 02, . . .) 
and substituting 

t' ^ t — a, t ^ t + a. 
p.2|) thereby takes the form 

r dz_^-2i{a,z),^ _^ ^ ^ [z~^])r(t - a - \z~^]) ^ 0. 

It is convenient to introduce the elementary Schur functions hn{t), n = 0, 1, . . . 
here. They are defined by an exponential generating function as 

00 

n=0 

One can thereby express the exponential factor e"^^^"'^) as 



a, 2 



^/i„(-2a)z". 



n=0 

Moreover, the product of two shifted tau functions can be expressed as 



T{t + a+ [z-i])r(t -a- [z^^]) = exp f + ^)^" ) ' 

00 

= ^/i„(D,)^-"e<"'^'>T(t)-T(t), 



n=0 



where 



n=l 
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(|2.2p can be thus eventually converted to the Hirota form 

OQ 

/i„(-2a)/i„+i(A)e<"^^*V(t) • r(t) = 0. (2.3) 

n=0 

Note that this procedure is reversible; one can trace it back and recover (|2.2p 
from (ESI). 

The last equation (j2.3|) is a generating functional form of an infinite number 
of Hirota equations, which are obtained by expanding (j2.3p in powers of a. For 
example, the terms linear in a„'s gives the special Hirota equations 

(i?iA.-2/i„+ipt)...)T(<)-T(t)=0, n = 3,4,.... (2.4) 

Note that those for n = 1 and n = 2 are trivial identities. The equation for 
n = 3 is exactly the lowest Hirota equation in (|2.ip . 

2.2 Fay-type identities 

The following identity, referred to as the Fay identity in the following, holds for 
all tau functions of the KP hierarchy: 

(Ai - A2)(A3 - A4)T(t + [Ar^] + [X^'Mt + [X^'] + [AI^]) 

- (Ai - A3)(A2 - A4)T(t + [Ar'] + [Xs'Mt + [A^-i] + [A^ i]) 
+ (Ai - A4)(A2 - A3)r(t + [Ar^] + [X^'Mt + [X^'] + [X^']) - 0. (2.5) 

Ai,A2,A3 and A4 are arbitrary parameters. This is a generalization of Fay's 
trisecant identities for the Riemann theta function and the prime form on an 
arbitrary Riemann surface. Sato and Sato [29j derived this identity and its gen- 
eralizations as a consequence of Pliicker relations on an underlying Grassmann 
manifold. One can derive it from the bilinear identity (j2.2p as well, thought we 
omit details. We shall show a similar procedure for the differential Fay identity 
below. 

Adler and van Moerbeke [27] obtained the differential Fay identity 

(A - /.)r(i + [A-i] + [/.-i])r(t) - (A - f,)r{t + [A-i])r(t + [/i-^]) 

+ (air)(t + [A-i])r(t + - (9ir)(t + [^l-'])Tit + [X-']) = (2.6) 

by a confluence procedure letting two of A's in the Fay identity tend to 0. We now 
derive this identity from the bilinear identity (|2.2p . Firstly, let us differentiate 
the bilinear equation by t[. This gives the equation 

Secondly, assuming that A and fi sit on the far side of the contour |z| = i? of 
the integral (i.e., |A| > R and > R), we specialize t' as 

t'=t+[A-i] + [/i-i]. 
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Then, by virtue of the identity 



n \ w / \ w 

n=l 



the exponential factor e^^* reduces to a rational function of z, A, as 

\ XJ \ fiJ (z-A)(^-At)' 
The foregoing equation thereby becomes 

dz 



27r« (z — A)(z — fi) 



\z-'] 



Now assuming that T{t) is a holomorphic function of < in a neighborhood of 
t — 0, we conclude that the integrand is a mcromorphic function of z in \z\ > R 
with poles at z — oo, A, /i. Therefore we can use residue calculus to calculate 
the contour integral as a sum of the residues at the poles. This yields (|2.6p . 

I|2.6p is not very convenient for some purposes. It is more suggestive to 
rewrite it as 

T(t+[A-i] + [/.-i]Mt) 1 r(t+[A^i]) 

r(* + [A-i])r(t + [/i-i]) A-M ^ ^T(t+[A.-i])- 

Moreover, one can use the differential operator 

oo „„ 
Z 



Diz) = Y.—9' 

L ^ TO 



Tl — 1 



to rewrite it further as 



exp I (e 



((e-(^) - l)(e-(^) - 1) log.) ^ 1 - (2.8) 
V / A — II 



As proven by Takasaki and Takebe jj, the differential Fay identity (|2.6p is 
equivalent to the KP hierarchy itself. A clue of the proof is a close relationship 
between the differential Fay identity and the auxiliary linear problem. 

2.3 Auxiliary linear equations for wave functions 

The auxiliary linear problem for the KP hierarchy [28l [29l [30l [31] consists of 
the linear equations 

(a„ - B„)* = 0, (2.9) 
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where i3„'s are differential operators of the form 

Bn - + bn,2dr^ + ■■■+ 6„,„. 

These operators (the Zakharov-Shabat operators) satisfy the Zakharov-Shabat 
equations 

[d,n-B,n,dn-B,,]=0, (2.10) 

equivalently, 

9„(B„) - d,n{Bn) + [B^, B„] = 0, 

where dn{Bm) means differentiating the coefheients of by i„. Moreover, 
there is a pseudo-differential operator (Lax operator) of the form 

oo 

such that the Zakharov-Shabat operators are expressed as 

Bn = (L")>o, (2.11) 

where ( )>o stands for the part of nonnegative powers of di. The Lax operator 
satisfies the Lax equations 

dn{L) ^[Bn,L]. (2.12) 

In analogy with quantum mechanics, solutions of the auxiliary linear equa- 
tions are called "wave functions" . The most fundamental wave function can be 
obtained from the tau function by the formula 

^(M) = l(^^e^(*^^) = f:^e«*-). 

T(t) T{t) 



As an immediate consequence of the bilinear equation (j2.2p for the tau function, 
this wave function and its dual 

^ ' ' ~ Tit) " Tit) ^ 

satisfy the bilinear equation 

^|^*(*',^)**(t,^) = 0. (2.13) 



We can derive the auxiliary linear equations (|2.9p from this bilinear equation 
as follows [33l [34] . Differentiating (|2.13|1 by t'^ yield the new bilinear equation 



/|^(5„*)(t',z)**(t,z) = 0. 
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In a similar manner, bilinear equations of the form 

/^(i?M')(t',z)vI/*(t,z) = 
J 27ri 

hold for for any higher order differential operator B = B{di) with respect to 
ti . (Actually, the bilinear equation (|2.13p should be interpreted as a generating 
functional form of all those equations.) Consequently, we obtain the equation 

(f^R{t,z)^*{t,z)^0 (2.14) 

for the linear combination 

R{t, z) = a„«'(t, z) - B^(t, z). 

We now choose i3 to be a differential operator Bn — Bn{di) of the aforemen- 
tioned form such that 

R(t, z) = 9„*(t, z) - B„ idi)^{t, z) = 0(z-i)e«(*'^), (2.15) 

namely, the Laurent expansion of i?(<, z)e~^*^*'^-' at z = oo is of order 0{z~^). 
Such an operator _B„(3i) does exist and is unique (see below for an explicit 
construction of this operator). On the other hand, there is a general formula 
that relates a contour integral of the form 

to the product of P{x, dx) and Q{x, dx)* (the formal adjoint of Q{x, dx) [551151] . 
With the aid of that formula, we can prove that any function i?(t, z) with these 
properties (|2.14p and (|2.15p has to vanish identically. Thus 5'(t, z) turns out to 
satisfy the auxiliary linear equations (|2.9p . 

Having obtained the auxiliary linear equations for ^(t, z), we can deduce 
the existence of the Lax operator as follows. A clue is the so called dressing (or 
Sato- Wilson) operator 

oo 

W = 1-h^WjV 

defined by the coefficients of the Laurent expansion 

oo 

«'(t,z)e-«(*^") = l + ^WjZ~^ 

of 1'(*,z)e-«(*'^) at z = oo. Note that 5'(/:, z), in turn, can be expressed as 

^'(<,z) = VFe«(*'^'. (2.16) 
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One can immediately see that the operators 



turn out to fulfill the condition p.lSp . Therefore if we define L as 

L = WdiW-^, 



the foregoing expression (|2.11|) of i3„'s is achieved. Moreover, by plugging (I2.16P 
into the auxiliary linear equations, one finds that W satisfies the the Sato equa- 
tions 

d^iW) = B^W - Wdl ^ -{Wd'lW^^)<oW, (2.17) 

where ( )<o stands for the part of negative powers of di. The Lax equations 
(|2.12p readily follow from the Sato equations. 

Note that the identification of the Lax and Zakharov-Shabat operators is in- 
dependent of the foregoing derivation of the auxiliary linear equations. Namely, 
once one can anyhow show that z) satisfies the auxiliary linear equations 
(|2.9p . one can thereby conclude the existence of an operator L which satisfies the 
Lax equations (|2.12p and by which i?„'s are expressed as (|2.1ip . This remark 
becomes technically important when we derive the auxiliary linear equations 
from the differential Fay identity. 



2.4 Differential Fay identity and auxiliary linear equations 

The foregoing derivation of the auxiliary linear equations (|2.9p from the bilinear 
equation (|2.13p is somewhat indirect. We now show that the differential Fay 
identity gives a more direct approach to the auxiliary linear equations. 

Let us anyway rewrite the differential Fay identity in the language of the 
wave function. The first step is to shift t in ^TT} as t ^ i - [A^^] - [/z^^]. ^TT^ 
thereby changes as 

r(f-[A-i]-[/x-i]Mf) 1 r{t-[\-^]) 

r(t-[A-i])r(t-[M-i]) ^ A - ^ ^ T[t - [^,-^])■ 

The next step is to multiply both hand sides of the equation by (A— /i)e^(*^^V(i — 
[/i~^])/T(t). This yields the equation 

(A-M)e^^*-^^e-^(^) ^^*~/^"'J^ 
r(t) 

= A + 9ilog^ l—^-^i-di\og 



r{t) r{t) 



X e 



{t - [^1- 



r{t) 
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Noting that the commutator of e ^'^^^ and e^'^*'''-' yields the factor 1 — /i/A, we 
can rewrite the last equation as 

Ae^^(^)*(<, /i) = (9i log A) - di)-^{t, fi). (2.18) 

Expanded in powers of A, this equation yields an infinite number of linear equa- 
tions of the form 

/i„(-(9t)«'(*,Ai) A/) (2.19) 
for n = 2, 3, . . ., where dt is defined as 

and VnS are the coefficients of the Laurent expansion of di log^(t, A), namely, 

oc 

di log4'(t,A) = A + ^i;„+iA-". 

n=l 

Linear equations (I2.19P were first discovered by Sato and Sato [29^ , but im- 
plications of these equations remained to be clarified for years. It was then 
pointed out by Takasaki and Takebe [J that (|2.19p are equivalent to the stan- 
dard auxiliary linear problem (|2.9p . As one can immediately see from the first 
few equations 

the left hand side of the n-th equation of ()2.19|) takes such a form as 

h,,{-dt)^(t, /i) = (^-ia,, + C„(9i, . . . , 9„_i)^ *(*, ^), 

where C„(9i,--- ,9„-i) is a differential operator that contains di, . . . ,dn~i- 
One can recursively eliminate the operators 82, ■ ■ ■ ,d„-i other than di using 
the preceding equations hm{—dt)'^{t,z) = Vrn^{t,z) for m < n. The outcome 
are linear equations of the usual evolutionary form 

^2^{t,^l) = {^f-2v2)^{t,^l), 

ds^it, //) = {df - 3v2di - 2diV2 - 3w3)*(<, /^), • • • • 

One can identify the operators on the right hand side to be the standard 
Zakharov-Shabat operators (|2.1ip by the same reasoning as we have used for 
deriving the auxiliary linear equations l|2.9p from the bilinear equation (|2.13p . 
Thus (|2.19p turns out to be equivalent to (|2.9p . In other words, (|2.18p is a 
generating functional form of these linear equations. 
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We are thus led to the conclusion that the differential Fay identity is actually 
the auxiliary linear problem in disguise. This fact lies in the heart of the proof 
[1] of the fact that the differential Fay identity is equivalent to the full system 
of the KP hierarchy. 

Let us note that v^s are rather well known quantities, namely, conserved 
densities of the KP hierarchy [35]. Since A) — A^'(t, A), one can rewrite 

the foregoing defining equation of u„'s as 

oo 

A) = (a + «„+iA-") *(t, A) 

CO 

= (l + ^?;„+iL-")*(«,A). 

ra=l 

This implies the operator identity 

oo 

Si =L + ^w„+iL-". (2.20) 

n=l 

Moreover, since 

log A) = A) + (e-^(^) - 1) logr, 
one has the explicit formula 

Vn = dih„{-dt)l0gT. (2.21) 

The differential Fay identity is particularly useful for studying the limit to 
the dispersionless KP hierarchy. This limit is achieved by a kind of "quasi- 
classical" limit. One can formulate this procedure in both the Hirota formalism 
and the Lax formalism. Let us now turn to this subject. 

2.5 Quasi-classical limit in Hirota formalism 

In the Hirota formalism, the quasi-classical limit is formulated as follows |36j . 
Firstly, we allow the tau function to depend on an extra small parameter H 
(an analogue of the Planck constant in quantum mechanics) as r = T(h,t). 
Secondly, we assume that the rescaled tau function Tuit) = T{h,h^^t) behaves 
as 

Tnit) = exp(ft-2F(t) + O(ft-i)) (2.22) 

in the limit of ft —s- 0. F(t) is a key function in the theory of dispersionless 
integrable systems, having several different names such as "dispersionless tau 
function" , "free energy" , "prepotential" , etc. that stem from its origin in prob- 
lems of mathematical physics [U El [10] . 
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By rescaling the tau function as above, the differential Fay identity 
takes the rescaled form 



exp 



A — /i 

Under the foregoing ansatz (|2.22p . both hand sides of this equation has a finite 
hmit as ?i — > 0. This yields the equation |4j 

^DWDi,)F ^ 1 _ d,{D{X)-D{„))F ^2.23) 

A — /i 

for the F function F — F{t). 

The last equation, known as the dispersionless Hirota equation, is a gener- 
ating functional form of an infinite number of equations, which are obtained 
by Laurent expansion of both hand sides in powers of A and /i. It should be 
noted that "the dispersionless Hirota equation" might be a misleading name, 
because (|2.23p is neither a direct limit of the Hirota equations (|2.ip themselves 
nor of their generating functional form (|2.2|) , but of the differential Fay identity. 
Therefore, it would be better to call (|2.23|) "the dispersionless differential Fay 
identity" . 

In this respect, Carroll and Kodama [12] observed an interesting fact: They 
pointed out that the special Hirota equations (|2.4p have a direct quasi-classical 
limit. To see this fact, it is convenient to start from the generating functional 
form 



of (|2.4p . It is not hard to see that this equation survives the quasi-classical limit 
and yields the equation 

oo 

eD(MfF ^^^^^-n^lg^g^p (2.24) 
n=l 

for the F function. Actually, the last equation can also be obtained by letting 
A ^ /i in (|2.23p . Expanded in powers of /i, this equation generates an infinite 
number of equations of the form 

dAF^K+,{Z,,Z2,...), n = l,2,..., (2.25) 

where 

d.dkF 



jk 



j+k—n 

Carroll and Kodama studied the algebraic structure of these genuine "disper- 
sionless Hirota equations" in detail. 

As the differential Fay identity is equivalent to the KP hierarchy itself, its 
dispersionless analogue (|2.23p turns out to be equivalent to (the Lax formalism 
of) the dispersionless KP hierarchy [141 HH |25] . This is by no means an obvious 
fact. We shall return to this issue later on. 
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2.6 Quasi-classical limit in Lax formalism 

In the Lax formalism, the procedure of quasi-classical limit starts from the WKB 
ansatz [371 ISH] 

z) = cxp(^h-^S{t, z) + 0(S°)) (2.26) 

of the rescaled wave function 'i'h{t, z) — ^'(fi., hr^t, z). It is easy to see that this 
WKB ansatz follows from the quasi-classical ansatz (|2.22p of the tau function. 
The associated S function is given by 

oc oo 

S{t, z) = at, z) - D{z)F = - E —dnF. 

n— 1 n— 1 

The auxiliary linear equations for the rescaled wave function takes the rescaled 
form 

hdn-^H{t,z) - Bn,n{hdi)-^h{t,z), 

where the coefficients of the operators 

BnA^di) = {hdiT + bKuAf^iT-'' + ■■■ + bn,n,u 

have a smooth limit as fi. 0. These equations may be thought of as analogues 
of the time-dependent Schrodinger equations in quantum mechanics. Conse- 
quently, the 5* function satisfies the Hamilton-Jacobi equations 

dnS{t,z) = Bn{diSit,z)) (2.27) 

with the Hamiltonians 

Bn{p) = lim Bh^p). 

These stuff resemble the setup of quasi-classical approximation in quantum 
mechanics. The Zakharov-Shabat operators -B;t.„ = Bfi^hdi) are now replaced 
by functions Bn = Bn{p) of a new variable p, which may be interpreted as the 
conjugate momentum of the coordinate x = ti. By this "quantum-classical cor- 
respondence" , commutators of differential operators in the Lax and Zakharov- 
Shabat equations of the KP hierarchy are replaced by Poisson brackets 

{F,G} = {dpF){diG) - {diF){dpG) 

of functions on a two-dimensional phase space. We are thus led to the Zakharov- 
Shabat equations 

dn{B,n) - d^Bn) + {B^M = (2.28) 
and the Lax equations 

a„(£)={S„,£} (2.29) 
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with respect to Poisson brackets [371 ISE SD] • The Lax operator is also replaced 
by a function of the form 

oo 

BnS are thereby expressed as 

Bn = (/:")>o, 

where ( )>o denotes the polynomial part of Laurent series of p. 

The Hamilton- Jacobi equations (|2.27p and the Zakharov-Shabat-Lax equa- 
tions l|2.28p - (|2.29p can be transferred from one side to the other by the func- 
tional relation 

oo _^ 

p = diS{t,z) = z-Y'- — diduF (2.30) 
^-^ n 

n—l 

connecting the spectral parameter z and the momentum p. The inverse p i-^ z 
of the map z i-^ p is exactly the Lax function C = C{p). One can see this fact 
from p.2ip and (|2.20p as follows. By rescaling the tau function, p.2ip becomes 

Vh,n = fldihn{-hdt)\0gTn,, 

hence 

lim vnn = didnF. 

h^o n 

Consequently, (|2.20p yields 

p = c{p)~-y^^^d,d^F 

n=l 

in the limit as fi. — ^ 0. Comparing this relation with (|2.30p . one finds that the 
inverse of the map z p — diS{t, z) is given by z = C{p). 

2.7 Dispersionless Hirota equation and Hamilton- Jacobi 
equations 

As already mentioned, the dispersionless Hirota equation p.23p is equivalent to 
the Lax formalism of the dispersionless KP hierarchy [TH [531 US] • Main part of 
the proof of this fact is to show that (|2.23p is a generating functional form of 
the Hamilton- Jacobi equations (|2.27p . This is a dispersionless analogue of the 
relation between the differential Fay identity and the auxiliary linear equations. 

The notion of "Faber polynomials" in complex analysis [25] plays a central 
role here. Let us demonstrate its power by assuming (|2.23p and deriving (|2.27p 
therefrom. We define p{z) as 

p{z) = diS{t,z) 
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and rewrite (|2.23p as 



,D(\)D(^,)F ^ p(A) -p(Af) 

X — IJL 



Taking the logarithm of both hand sides yields 
By virtue of the identity 



we can further rewrite the last equation as 

n—l \ m—1 / 

^-^ n 

n— 1 



(2.31) 



The proble is to solve this equation for dnS{^). To this end, we introduce the 
Farber polynomials <i>„(p), n = 1, 2, . . ., of the map z p{z), which are uniquely 
defined by the generating functional relation 

log — = - > 

z ^ — ' n 

n—l 

More explicitly, as we shall show below, ^nip) are given by the polynomial part 
of the inverse map p t-^ z = z{p) as 

$„(p) = (z(p)«)^^. (2.32) 

In the present setting, z(p) is nothing but C{p), hence 

$„(p) =S„(p). 

We can now solve (|2.3ip for 9„5(^) as 

dnSifi) - $„(p(a*)) = Bnipiix)). (2.33) 

The last equations are exactly the Hamilton- Jacobi equations (j2.27p . 

Let us show that (|2.32p indeed holds. To this end, we differentiate the 
defining equation of <i>„(g) by z. This yields the identity 

U \ °° 
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where p'{z) denotes the z-derivative dzp{z). We can extract $„(g) by a contour 
integral as 



where the contour of integral is a sufRciently large circle \z\ = R such that 
|g| < R. Since the map z i-^ p{z) is invertible in a neighborhood of z = oo, 
we choose the contour inside that neighborhood and change the variable from 
z = z{p) to p as 



The contour of the transformed integral is a closed curve C that encircles p — oo 
but not q. This integral is nothing but the polynomial part of z{q)"'. 

3 Toda hierarchy and difference Fay identity 
3.1 Bilinear equations for tau function 

The Toda hierarchy [41] has a discrete variable s £ Z (lattice coordinate) and 
two sets of continuous time variables t — {ti,t2, . . .), t ~ {ii,t2, ■ ■ ■)■ The tau 
function t — t(s, t, t) satisfies an infinite number of Hirota equations 



where Dn and £)„ are the Hirota bilinear operators associated with 9„ — d/dtn 
and Bn — d/dtn- A generating functional form of these Hirota equations is given 
by the bilinear equation 



Whereas the contour of integral on the left hand side is a sufficiently large circle 
\z\ = R, that of the right hand side is a sufficiently small circle \z\ = R~^ around 
the origin (or one may understand both integrals as a purely algebraic operator 
extracting the coefficient of z^^). 

It is easy to from this bilinear equation that the Toda hierarchy contains the 
KP and modified KP hierarchies [30] as subsystems. For example, if t' and s 
are specialized asi' = i and s' = s, (|3.2p reduces to the bilinear equation 





DiDit{s, t, t) ■ r(s, t, t) + 2t(s + l,t, t)T{s - 1, t, t) = 0, 
(£12 + DI)t{s + l,t~t)- r(s, t, t) = 0, 
{D2 + Dl)T[s,t~t)-T[s + l,t~t)^Q, 



(3.1) 




(3.2) 
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which is substantially the bilinear equation for the KP hierarchy. Thus the tau 
function of the Toda hierarchy, viewed as a function of t, is a tau function of 
the KP hierarchy as well. The last equation can be generalized to the bilinear 

equations 

/^e«(*'^*'"V(s',t'-[z-i],<)T(s,t+[z-i],t) = for s'>s, (3.3) 
J 27rj 

which show that the tau function may be thought of as the tau function of the 
modified KP hierarchy with respect to t and s. The same interpretation holds 
true for i and s. 

As in the case of the KP hierarchy, one can rewrite this bilinear equation 
into the Hirota form 

OG 

J2 /i„(-2a)/i„+,._,+i(A)e<"'^'>+<^'^*>T(s, t, t) ■ t(s', t, t) 

oo 

= ^/i„(-2a)/i„_,-+,_i(^i)e<"'^*>+<"'^J>r(s- -rGs' + l,*,*), (3.4) 

where a = (01,02,...) and a = (01,02,...) are new variables, and and 
(a, Di) are the counterparts of Dt and (a, Dt) for i, namely, 

D-t = (z)i, . . . , . . (a, A) = ^a„Z)„. 

n=l 

Expanded in powers of a and a, (|3.4p generates an infinite number of Hirota 
equations. 

3.2 Difference Fay identities 

We can derive three Fay- type identities [111132] from p.2p by specializing s', t', t' 
as follows. 

(i) s' ^ s + \,t' = t+ [A-i] + [^-1], i = t. 

(ii) s' = s - 3, t' = t, = t + [A] + 

(iii) s' = s, t' = t+ [A-i]. t' =t+[ii]. 

As we have shown in the case of the KP hierarchy, the exponential factors in 
the bilinear identity thereby become rational functions, and calculation of the 
contour integrals reduces to residue calculus. 

Note that the role oi x = ti in the KP hierarchy is now played by s. Unlike 
that case, however, we now shift s rather than differentiate the bilinear equa- 
tion. For this reason, we call the following Fay-type identities "difference Fay 
identities" . 
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(i) The specialized bilinear equation becomes 

/dz - 
— zt(s + 2, t, t - [0])r(s - 1, t, * + [^]). 

Residue calculus yields the first difference Fay identity 

T{s + l,t+[X-'] + [ii-'])T{s,t,t) 

^ T{s + l,t+[iJ,-%i)T{s,t + [X-^],i) 



\ — ji 



+ i-^t{s + l,t+ [\'\t)T{s, t + [ix-\t) = 0. (3.5) 

A fJi 

(ii) The specialized bilinear equation becomes 

^^" (.-A)(.-M) "^^ - 2, t, t + [A] + M - [.])r(. - 1, t, * + [.]). 
Residue calculus yields the second difference Fay identity 

^T{s-2,t~t+[\] + [lj\)T{s-l,t,t) 

+ x{X-^f ^^ -2,t,i+ [fi])T{s -l,t,i+ [A]) 

^ T{s-2,t,i+[X])T{s-l,t,t+[ij]) = 0. (3.6) 



(iii) The specialized bilinear equation becomes 

= / ^t(s + 1, t + [A-i], * + [m] - [z])t(s - 1, t, t + [z]). 
J Zm z — n 

Residue calculus yields the third difference Fay identity: 

XT{s,t + [X-%i+[pi])T{s,t,i) 

- Xt{s, t,i + [X-%i)T{s, t,i+ M) 

+ ht{s + l,t + [X-%i)T{s -l,t,i + [m]) = 0. (3.7) 



19 



These difference Fay identities can be cast into different forms. The following 
are counterparts of (|2.7p : 



t{s, t, t+[X] + [m])t(s, t, t) _ t(s, t, t) 

,1,1)'' 

, (3.9) 



-{s, t, t + [X])t{s, t, i + [^i]) (A - ^)t(s + 1, t, t) 

^lT{s+l,t,t+[X]) XT{s + l.t,t+[^i]) 



T{s,t,t+[X]) T{s,t,t+[^i]) 
T{s,t+[X-\t+[^i])T{s,t,t) 

T[s,t + [x-^]:t)T{s,t~t + [^Ji]) 

_ ^ /i t(s + 1, t + [X-\1)t{s ~l,t~t+ 

X T{s,t + [x-^,t)T{s,t,t + [^l]) ■ ^ ■ > 



Similarly, the following equations are counterparts of 
exp((e^(^)-l)(e^(^)-l)logr) 

= ^-^(Aexp((e^(^) - l)(e-^= - l)logT) 

- Mexp((e^('') - l)(e-^= - l)logr)), (3.11) 

exp((e^(^) - l)(e^(^) - l)logr) 

= ^(-/iexp((e^(^) - l)(e^= - l)logr) 



+ 



Aexp((e^('') - l)(e^-' - l)logT)), (3.12) 



exp 



,(^(e^(A) _ i)(e-S(f) - l)logr) 
= 1 - ^exp((e-°(^) - l)(e^= - 1) logr + (e-^'^^ - l)(e"^= - l)logr 

-(e^= -l)(e-^= -l)logr), (3.13) 

where 



oo _^ oo ^ 

D(z) = V— 9„, 5(z) = V— 9„. 

JL / JL / 

n—1 n—1 



20 



3.3 Auxiliary linear equations for wave functions 

We now introduce the wave functions 

T(S,t,t) 
T{s,t,t) 



and the duals 



T{s,t,t} 
T[s,t,t) 



The bilinear equation p.2p for the tau function becomes the bilinear equation 

/ — ^-rs', t', i, z)^*(s, t, lz)= / — *(s', t', t\ z)^*{s, t, i, z) (3.14) 
/ 27r? / 27ri 

for these wave functions. From this bilinear equation, one can derive the auxil- 
iary linear equations 

(5„-B„)$ = 0, (4-B„)$ = (3.15) 

for $ — '^{s, t, i, z), ^{s, t, i, z). Bn and i?„ are difference operators of the form 



and satisfy the Zakharov-Shabat equations 

[dm — Bn, dn — Bn] = 0, [5,„ — i?„, 5„ — i?„] = 0, 
[dm — Bm, On — Bn] = 0. 

We can construct two Lax operators of the form 



(3.16) 



j=o j=i 



such that the Zakharov-Shabat operators are expressed as 

Bn = (i")>o, Bn = (i")<0, 

where ( )>o and ( )<o are the projection onto nonnegative and negative pow- 
ers of , respectively. The Lax operators satisfy the Lax equations 

dn{L)^[Bn,Ll dn{L) =[i?„,L], 

dn{L) = [Bn,L], dn{L) = [5„,L]. ^ ' ' 
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3.4 Difference Fay identities and auxiliary linear equa- 
tions 

We now rewrite the difference Fay identities (|3.8|) - (|3.10p in the language of the 
wave functions. This leads to a generating functional expression of the auxiliary 
linear problem. 

Let us first consider (|3.8p . We shift the time variables a,st-^ t—[X^^] — [fJ,~^] 
and muhiply both hand sides by {\- ^i)e^^*'^'^ /t(s + 1, t, t)T{s,t-[\-^],i). This 
yields the equation 



T(5,f-[/i-i],t) 
r{s,t,t) 



^ ^^ rjs + 1, t - [X-^],t)/T{s + 1, t, t) t{s, t - [ti-\t) 

T{s,t-[X-^],t)/T{s,t,i) ^ Tis,t,t) 

r{s,t,t) 

Noting that the commutator of e^'^^'^^ and e^*^*'''-* yields the factor 1 — /i/A, we 
can rewrite the last equation as 

Ae-^(^)vl.(s, t, t, ^i) = ^if±l^^vl;(s, t, Ai) - + 1, t, t, (3.18) 
^(s,t,t, A) 

Let us now consider (|3.9p . We now shift i and s as t ~> t — [A] — [/i] and s 
s + 2, and multiply both hand side by A(A — ^)^'*e^^*'^ /T{s,t,t)T{s,t,i— [A]). 
This yields the equation 

T{s,t,t) 

^ 4(t,,,-i)r(s + MJMM:^ _ _^ T(g + l,f,f-[A])/r(.,t,t) ^^ 

T(s,t,t) T(s,t,t-[A])/T(s-l,i,t) 

X .-ig?(t,M)l(vMMMD 
r(s-l,t,t)' 

which can be cast into the form 

e"^(^)'^(s, t, t, ^i) = ^-(5, i, t, - - - 1, t, t, ^). (3.19) 

w(s — 1, t, A) 

In the same way, we can derive two more linear equations for the wave 
functions from (|3.7I1 . Firstly, if we shift the time variables ast^t — [X~^] and 
and multiply both hand side by X''e^^*''^yT{s,t,i)T{s,t,i — [/x]), we 
obtain the equation 

e-^^^^^{s, t, I A) = ^{s, t, t, A) - J^(^^ M^^) vi/(, _ 1, t, I A). 

W(s — 1, t, t, II) 
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It will be better to exchange A and /i as 

e-^(^)^'(s,t,<,M) = ^{s,t,i,fi) - J^''*'l'^\^ 'f{s-l,t,t,fi), (3.20) 

because the outcome takes the same form as p.l9p . Secondly, if we shift t,i 
and s as t —> t — [A~^], i —> i — [^], s —> s + 1 and multiply both hand side by 
X^'^e^'-^''^ /t{s + l,t,i)T{s,t — [A^-^],t), we are led to the equation 

Ae-^(^)'J'(s, t, I f,) = ^lif±llSlL^*(s, t, I f,} *(s + 1, t, t, y,). (3.21) 
^(s,t,t, A) 

p.lSp - (|3.21|) show that the linear equations 

Ae--(^)$ = ^^(f±l:SMa>-e^=$, 

^'(s,t,t,A) 

e-DW^ = $ - JfijM^^-a,^ 
*(s- l,t,<,A) 

are satisfied by both $ = ^'(s, t, t, /^) and $ = ^'(s, t, y). Expanded in powers 
of A, they generate an infinite number of linear equations of the form 

K{-dt)^^v^^, /i„(-9i)$ = -5„e-^=$, (3.22) 

where dt is the same notation as defined for the KP hierarchy, di is similarly 
defined as 

d2 d„ 



K fa "2 On 

V 2 n 



and Vj = Vj(s,t,t) and Vj = Vj{s,t,t) are the Laurent coefficients of the ^'- 
quotients in the foregoing linear equations, namely, 

^(s + l,t,t,A) _ , , ^ ,i_„ §(s,f,f,A) ^„ 

Note that the lowest (n = 1) equations of p.22p read 

9^$ = (e^= - wi)$, ai$ = Uie-^=$. (3.23) 

As we have illustrated in the case of the KP hierarchy, one can derive the 
usual evolutionary form (j3.15p of auxiliary linear equations from these equations 
p.22p and (|3.23p . On the basis of this fact, Teo |32] proved that the difference 
Fay identities are equivalent to the full system of the Toda hierarchy. 

3.5 Dispersionless Hirota equations and Hamilton- Jacobi 
equations 

The quasi-classical ansatz p.22p for the tau function of the KP hierarchy can 
be readily generalized to the Toda hierarchy [l^- Namely, we allow the tau 
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function to depend on h as t = T{h, s,t,t), and assume that the rescaled tau 
function Tri{s, t, i) — T{h, h^^s, hr^t, hr^i) behave as 



(s, t, t) = exp(h-^F{s, t, t) + 0(ft~^)) (3.24) 



in the hmit oi h ^ Q. Note that the rescahng changes the lattice spacing from 
1 to /i, eventually tending to in the quasi-classical limit. The discrete variable 
s in r(/i, s, t, t) is thereby replaced by a continuous variable in i^(s, t, t). 

Under this quasi-classical ansatz, we can derive the following three disper- 
sionless Hirota equations [20l[2Tl[T4] for the F function F = F{s,t,i) from the 
difference Fay identities: 

D{\)D{p)F _ ^6 ~ Me 

\ ' 
A — fj, 



gS(A)5(p)F ^ -/ie-^-^--- +Ae-^^'^--- (3.25) 



A — /i ' 

^D{X)D{t,)F A^g(D(A)-5(p)+as)a.F^ 



A 

It will be better to rewrite the second equation as 

A-i - /i-i 

by which the symmetry between t and t becomes manifest. Moreover, defining 
the S functions as 

S{z) ^ i{t,z) ^ s\ogz - D{z)F, 
S{z) = ^{i; z-^) + s log z + dsF- D{z)F, 

we can rewrite p.25p as 

A — /i ' 

^^g-a,s(A)_g-a.s(p) (3 26) 



gD(A)D(p)F 



5(A)F)(p)F _ d'^F 

^D(\)D(^t)F ^ ^ _ g-a,S(A)ga,S(M)^ 

The S functions are actually the phase functions in the WKB ansatz 
«'s(s, t, t, z) = eyip{h-^S{z) + 0(?i°) 



^'fi.(s,t,t,z) = cxp[h-^S{z) + 0(fi°)) 



(3.27) 



for the rescaled wave functions. The auxiliary linear equations thereby yield the 
Hamilton-Jacobi equations 

d^S = 6„(e^=^), Br^S = 6„(e^=^) (3.28) 
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for S = S{z),S{z). Here S„(P) and Bn{P) are classical counterparts of i?„ = 
B„(e^=) andB„ = S„(e^=): 

Bn{P) = lim Sr,„ (F) , ^„ (P) - lim Pft.„ (F) . 

The dispersionless Hirota equations in the form of (|3.26p turn out to be equiv- 
alent to these Hamilton- Jacobi equations, hence to the dispersionless Toda hi- 
erarchy itself [H [H US] . 

4 BKP hierarchy and differential Fay identities 

4.1 Bilinear equations for tau functions 

In this section, we mostly consider the two-component version [441 [45] of the 
BKP hierarchy [1^ . This system contains two copies of the one-component BKP 
hierarchy as subsystems. The situation is thus somewhat similar to the Toda 
hierarchy, which may be thought of as the charged two-component KP hierarchy 
as well |41[ . This is a reason why we are interested in the two-component case 
rather than the one-component BKP hierarchy. 

The two-component BKP hierarchy has two sets of time variables t = 
(ti, is, . . .) and i — (fi, ta, . . .) indexed by odd positive integers. The tau func- 
tion T = r(t, t) satisfies the bilinear equation 

/ ^S''-'''\{t' - 2[z-\t')T{t + 2[z-\t) 

J ZTTIZ 

= ^^e«(*'-*-^)r(<',r-2[z-i])r(t,t + 2[z-i]). (4.1) 

The integrals on both hand sides are understood to be a contour integral along a 
sufficiently large circle |z| = i? (or just a formal algebraic operator). The other 
notations are as follows: 

C(t, Z)=Y. t2n+lz'-+\ H ^ (a, ^, . . . , . . .) . 

n— 

If t is specialized as i' — t, the bilinear equation ([4.ip reduces to 




e«(*'-*^^V(t' - 2[z-^],t)T{t + 2[z-^],t) = T{t'~t)T{t~t). (4.2) 



Similarly, specializing t' as t' — t yields 

/ ^g«(t'-t,^)^(^^ I' _ 2[z-i])r(«, t + 2[z-i]) = r(t, <>(<, t). (4.3) 
J 2'Kiz 

These equations coincide with the bilinear equation of the one-component BKP 
hierarchy. Thus the tau function of the two-component BKP hierarchy may 
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be thought of as the tau function of the one-component BKP hierarchy with 
respect to both t and i. 

(|4.ip can be converted to the Hirota form 



oo 

= J2 /i4-2a)/i42A)e^"''''>+<"'''*>T(t, t) ■ T(t, t), (4.4) 

where a = (01,03,...) and a — (01,03,...) are new variables, Dt,Di, and 
(a, I?t), (a, Di) are operators defined as 

Dt = {Di, . . . , . . , (a, Dt) ^ 02„+li:'2n+l, 



n=0 

00 



n— 

and hn{t) are defined by the generating function 

00 00 



^) = exp(^t2„+i^2n+l ^J2hr,it)z^ 



n=0 n=0 

4.2 Differential Fay identities 

We can derive four differential Fay identities [3T from (|4.ip as follows. 

(i) Differentiate the bilinear equation by t[ and specialize t' ,i' as t' = t + 
2[A-i] + 2[/i"i],t' = t. 

(ii) Differentiate the bilinear equation by t[ and specialize t' ,t' as t' = t, 
r = t + 2[A-i] + 2[/i-i]. 

(iii) Differentiate the bilinear equation by t[ and specialize t, t' as i' = t + 
2[A-i], r = i + 2[/i-i]. 

(iv) Differentiate the bilinear equation by t[ and specialize t, i' as t' = t + 
2[A-i], t' = < + 2[/i-i]. 

Note that the exponential factors in the bilinear equations become rational 
functions by virtue of the identity 



^ 2 /z\2n+i_ l-z/w 
-H; 2n+ 1 \w) ^ °^ f + z/w' 



n=0 



For example, in the case of (i), 

(z - A)(z - ^i)' 

We can thus calculate the contour integrals by residue calculus. 
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(i) After differentiating by t[ and specializing t' , the bilinear equation becomes 

{d,T){t + 2[A-i] + 2[^^-\t))T{t + 2[z-'],t) 

{diT){t + 2[A-i] + 2[/x-i], t - 2[0-i])r(t, t + 2[z-']). 



+ I 

dz 



2niz 

By residue calculus, we obtain the first differential Fay identity 



T{t,t) 

T{t + 2[X-']+2[i^-%t)T{t,i) 
T(t + 2[A-i],t)T(t + 2[/i-i],t) 

^ + ^^-M-a.logt + y:!!'Sl- (4.5) 



A-mV ^^r(< + 2[/x-i],t) 

(ii) This case is substantially the same as (i), only the role of t and i being 
exchanged. We thus obtain the second differential Fay identity 

T{t,t) 

T{t,i + 2[X-^]+2[^x-^])T{t,t) 



r(t,t + 2[A-i])r(t,t + 2[M-i]) 

= A±if A - M- a. log t' ' + (4.6) 

(iii) After differentiation and specialization, the bilinear equation becomes 

+ {dir){t + 2[A-^] - 2[z-%i + 2[M-']))T(i + 2[z-^],i) 
= f ^ + 2[A-^],t + 2[/x-i] - 2[z-'])T{t,t + 2[z-']). 

By residue calculus, wc obtain the third differential Fay identity 

A-a.logli^+2[^"]'* + ^[^"^J)^ 



Tit,t) 

T{t + 2[X-^],i + 2[ii-^])T{t,t) 
r(t + 2[A-i],t)T(t,t + 2[M-i]) 



= A-a.logl|tM. (4.7) 
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(iv) This case is parallel to (iii) , and leads to the fourth differential Fay identity 

r"^^^°' — ^) — J" 

Tit + 2[X-^it + 2[^^-^])T{t,i) 
T(t + 2[A-i],t)T(t,t + 2[^-i]) 

We can use the BKP version 

°° ^-2n-l _ °o ^-2n-l _ 

n— n— 

of the D{z) operator to rewrite the foregoing differential Fay identities as follows: 



(^X + fi- ai(e2^(^)+2^W - l)logr) X 

xexp((e2^(^)-l)(e2^('')-l)logr 



= ^ (a - M - 5i(e2^(^) - e2^(^)) logr) , (4.9) 
A — /J, \ / 



(x + ^i- di{e^'^W+2DM _ 1) logr) x 
X exp 



= A±if (a - M - 5i(e25(^) - e2^('^)) logr) , (4.10) 
A — fj. \ / 



(A-ai(e2^(^)+2^('')-l)logT)x 

X exp((e2^(^) - l)(e2^('^) - l)logT) 

= A-5i(e2^(^)-e2^(''))logT, (4.11) 

(/x-ai(e2^(^)+2^('^)-l)logT)x 

X exp((e2^(^) - l)(e2^('') - l)logT) 

= M - ai(e2^W - e2^(^)) logr. (4.12) 
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4.3 Auxiliary linear equations 

We now introduce the two wave functions 

T{t,t) 

^t,t,z) = lM_^^Si.^). 

T[t,t) 

Their duals are given by '^'{t, t, — z) and ^{t, i, —z). The bihnear equation (|4.ip 
for the tau function yields the bilinear equation 

/ ^^{f, i', z)^it, t, -z) = (f ^^{f, t\ z)^it, I -z) (4.13) 
J Zmz J Zmz 

for these wave functions. From this bilinear equations, one can derive an infinite 
number of auxiliary linear equations of the form 

{d2n+l - B2n+l)^ = 0, (Szn+l " ^2„+l)$ = (4.14) 

for n = 0, 1, . . . and 

{didi-u)^^0 (4.15) 

that hold for both $ = ^{t,i,z) and $ = z). -B2n+i = -B2n+i(9i) and 

B2n+i — B2n+i{di) are differential operators with respect to ti and ti of the 
form 

B2n+1 = 9^"^^ + 62ri+l,2l9i" """ + ••• + &2ri+l,2nl9l, 
Bln+1 = B^"^^ + 62n+l,2^^" "'" + ••■ + &2n+l,2nl9l- 

Note that they do not have a 0-th order term. This is a characteristic of the 
operators that show up in the auxiliary linear problem of the one-component 
BKP hierarchy |46] . The first two sets of equations (j4.14|) are thus nothing but 
the auxiliary linear problem of the underlying one-component BKP hierarchies. 
The third equation (|4.15p may be thought of as the two-dimensional "integrable 
Schrodinger equation" associated with the Novikov-Veselov equation [48]. The 
potential u is given by 

u = -29i9ilogT. (4.16) 

In this respect, the two-component BKP hierarchy is also referred to as the 
Novikov-Veselov hierarchy, ti and ti play the role of spatial variables therein. 
Auxiliary linear equations of this type play a central role in Krichever's recent 
work [ini on a Schottky-type problem on Prym varieties [50] . 

The Lax formalism of this hierarchy exhibits new features because of the 
presence of two spatial dimensions. For example, the Zakharov-Shabat operators 
of the same type satisfy the usual zero-curvature equations 

[d2m+l — B2m+1, d2n+l — B2n+l] = 0, 
[d2m+l — B2m+1, d2n+l — -S2n+l] = 0, 
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but the equation for the pair of B2m+i and B2n+i has an extra term on the 
right hand side as 

[d2m+l - B2m+1, 82,1+1 ^ -82,1+1] = ^^mn (9l , 9l ) (^i 9i - u) , (4.18) 

where Dmn(di,di) is a differential operator with respect to both ti and ti 



4.4 Differential Fay identities and auxiliary linear equa- 
tions 

As we have demonstrated for the cases of the KP and Toda hierarchies, the 
differential Fay identities in the present case, too, give a generating functional 
expression of the auxiliary linear problem. The situation, however, turns out to 
be more complicated than the previous two cases. 

Let us first consider (|4.5p and (|4.6p . As regards (|4.5p . we first shift t as 
t ^ t — 2[A^^] — 2[/i^^]. The equation thereby changes to 

A + , + log J ^(i^^/" 

T(t-2[A-i]-2[/^-i],t)r(*,t) 



T(t-2[A-i],<)T(t-2[^-i],t) 



^ + '^^-M + 5ilog^(*-2[^'^]'*) 



A-/i V " "T(<-2[M-i],t) 

We then multiply both hand sides by e^(*'^'T(< — 2[/i^^],t)/r(/;,t). After some 
algebra, we obtain the equation 

{di log ^(t, t. A) + ai)e-2^(^'*(t, t, 

= (9ilog*(t,t, A) - (4.19) 

In much the same way, shifting i as, i ^ i — 2[A~^] — 2[/i~^] in (|4.6|) and 
multiplying it by e^'*'^V(/:,t — 2[^^^])/T(t, t), we obtain the equation 

(5i \og^{t, t, A) + ai)e-2S(^)5-(t, t, ^) 

= (ailog^'(<,<,A)-ai)*(t,t,M). (4.20) 



We can rewrite ()4.7p and (|4.8p in a similar way. As regards (|4.7p . we shift 
both t and i as t ^ t — 2[A~^] and i ^ t — 2[/i^^], and multiply both hand 
sides by e^(*''''T(t, t - 2[fi-^])/T{t,i). This leads to the equation 

(5i log ^(t, t, A) + ai)e-2^(^)*(t, t, Ai) 

= (ailog*(<,<,A)-ai)*(t,t,Ai). (4.21) 



Note that this equation has the same form as (|4.5[) . Anticipating a similar result, 



we now consider (14.81). If we shift t and < as f ^ f — 2fA ^1 and — > t — 2^" ^ 
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and multiply this equation by e^'-*''''V(t — 2[A ^]^t) /rit^t), the outcome is the 
equation 

(91 log ^{t, t, ^) + ^l)e~^^^^'^■^{t, i, A) = (^i log ^(t, i, fi) - d^)^(t, t, A), 
which is almost what wc want. By exchanging A and /i, we obtain the equation 
(9i log ^{t, t, A) + 9i)e-2C(A)^(^^ I ^) 

= (ailog§(t,t,A)-9i)*(t,t,Ai). (4.22) 

As anticipated, this equation has the same form as (|4.20p . 

These results show that '^{t^ i, fi) and ^{t, i, /i) satisfy the same linear equa- 
tions 

(5i log*(t, t, A) + ai)e-2^(^)$ ^ {di log*(<, i, A) - 

{di log*(t, t, ^) + ai)e-2S(^)$ = (Qi log*(t, /i) - di)^. 

Expanded in powers of A, they give rise to an infinite number of equations of 
the form 

(^hn+i{-2dt) + «„+i/i„_„(-2at) +5i/i„(-2at))$ = 0, 

rn—l 

for n — 2,3, . . ., where dt and di denote the set of operators 

a /a ^3 92,1+1 \ / 5 ^3 d2n+l \ 

5t = (5i,-,...,^,...j, a,= (9i,-,...^,...j, 

and VnS and i;n's are defined by generating functions as 

oc 

ailog^(t,t,A) = A + ;^«„+iA-", 

oo 

5i§(t,t,A) =A + ^^;„+iA-". 
For example, the equations for n = 2 in (|4.23p take the form 

(-^^3 + + 4{df \ogT)di) $ = 

or, equivalently, 

(93"93-6(a2logT)ai)$-0, 
(93 - df - 6{df logT)9i)$ ^ 0, 
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which coincide with the lowest equations of (|4.14p . The equations for n = 2k in 
(|4.23p take such a form as 



{-^^^2k+l - (5^(91, ... , d2k^l)^ = 0, 

where Cm{di, ■ ■ ■ , d2k-i) and Cm{di, . . . , d2k-i) are differential operators con- 
taining 9i, . . . , d2k-i and di, . . . , d2k-i, respectively, but no 0-th order term. As 
explained in the case of the KP hierarchy, one can convert these equations to 
the usual evolutionary form (I4.14p . 

Still missing is the two-dimensional Schrodinger equation (|4.15p . To derive 
this equation, we convert (|4.7p and (14. 8p to yet another form as follows. As 
regards (|4.7p . we first shift t and i as t ^ t — 2[A~^] and i ^ i — 2[/i~^], them 
multiply it by e^(*'-^V(t - 2[\-^],i)/T{t,t). This yields the equation 

(9i + di \og^{t, t, ^))e-2^(^)*(«, t, A) = {di - di log t, i, A). 

Exchanging A and /i, we eventually obtain the equation 

(5i -I- di log t, A))e-2^(^)*(i, t, fi) 

= (ai-5ilog^(t,t,A))*(t,t,Ai). (4.24) 

Similarly, shifting t and i as t ^ t — 2[A^^] and f — > f — 2[/i] in (14. 8p and 
multiply it by e^(*''^V(t,t — 2[/i])/T(t, t), we obtain the equation 

(5i + B, log t, A))e-2^(^)*(t, t, fi) 

= {^l-^^log^{t,lX))^{t,t,^^). (4.25) 

Expanded in powers of A, these equations give rise to an infinite number of 
linear equations. Among them, the equations from the A^^-terms give exactly 
the two-dimensional Schrodinger equation (|4.15p : 

(9i 91+2(^1 91 log T))^'(t,<,Ai) = 0, 
{didi+2{didi logr))^(t,t,Ai) = 0. 

Thus the usual auxiliary linear equations (|4.14p and (|4.15p can be derived 
from the differential Fay identities, ft should be noted that we have not fully 
used the generating functional equations (|4.19p - (|4.22p . (|4.24p and (|4.25p . Namely, 
only even power terms of A in the first four equations and the A^^ -terms in the 
last two equations are enough to recover all relevant auxiliary linear equations. 
What meaning, then, do other terms have? The fact is that the other terms 
give no new equations. For example, the equations for n = 3 in (|4.23p can be 
derived from the the corresponding equation for n = 2 by differentiating by ti 
and ti, respectively. 
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4.5 Dispersionless Hirota equations and Hamilton-Jacobi 
equations 

The procedure of quasi-classical limit is a straightforward generalization of the 
case of the one-component BKP hierarchy [ST] , 

As regards the Hirota formalism, we assume the quasi-classical ansatz 

Tr,(t, t) = exp(^h-^F{t, t) + Oih-^)) (4.26) 

for the rescaled tau function Tti{t,i) = T{h,h~^t,h~H). By this ansatz, we 
can derive the following dispersionless Hirota equations from the foregoing four 
differential Fay identities [17] : 

(a + /i - 2di{D{X) + D{fi))F^ ^iDWDMF 

= ^{^-f^-mD{X)-D{f,))F), (4.27) 



A + /i 



(^(^) - D{^i))F) , (4.28) 
(a - 2di{D{X) + 5(^))F)e4^(^)^(^)^ = A - 2di{D{X) - D[ii))F, (4.29) 

(ji - 2di{D{X) + 5(Ai))F) e4^(^)S(A')J^ = ^ _ 2di{D{n) - D{X))F. (4.30) 

By defining the S functions S{z) = S{t, t, z) and S{z) — S{t, t, z) as 

S{z) = C(t, z) - 2D{z)F, S{z) ^ i(t, z) - 2D{z)F, 
((07)) - can be rewritten as 

(ai5(A) + ai5(M))e4^(^)^(^)^ = ^(a,5(A) - d^Sifi)), (4.31) 

A — /i 

id,S{X) + 5i5(M)e4^(^)^(^)^ = ^{d,S{X) - d,S{p)), (4.32) 

A — /i 

(5i5(A) + 5i5(M))e4^(^)^('')^ = ai5(A) - diSifi), (4.33) 

(ai5(Ai) -)- 5i5(A))e4^(^)^(^)^ ^ ^lS{^l) - diSiX). (4.34) 
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Actually, (|4.3ip was first obtained by Bogdanov and Konopelchenko [35] by the 
9-dressing method, and rederived, along with the other equations, by Takasaki 
[17] by the method presented here. Later on, Chen and Tu [S3| studied these 
equations by the kernel method of Carroll and Kodama [H] . 

Quasi-classical limit is also achieved in the Lax formalism. The ansatz for 
the wave functions again takes the WKB form 



(4.35) 



^n{t, t, z) = eicp[h-^S{z) + 0{h°) 
^n{t, i, z) = cii^{h-^S{z) + 0{h")^ 



for the rescaled wave functions t, z) = '^{h, h^^t, h^H, z) and ^h{t, i, z) = 
^{h,h~^t,h~H,z). The foregoing S functions S — S{z),S{z) show up here as 
the phase functions of the WKB approximation, and satisfy the Hamilton- Jacobi 
equations 

d2n+lS — B2n+l{dlS), d2n+lS — B2n+l{dlS), , . 

OlOib — u 



with the classical Hamiltonians 



I32n+lip) = lim Bti,2n+lip), B2n+lip) = hm Bti,2n+lip) 

obtained from the (rescaled) operators Bf^^2n+i = Bnan+iif^di) and Bnan+i = 
Bh,2n+i{fidi)- Hamilton- Jacobi equations of this type were also studied by 
Konopelchenko and Moro [Sl] in the context of nonlinear optics. 

Though the structure of the dispersionless Hirota equations look consider- 
ably different from those of the dispersionless KP and Toda hierarchies, one can 
anyhow show that they are, in fact, equivalent to the Hamilton-Jacobi equations 



5 DKP hierarchy and differential Fay identities 
5.1 Bilinear equations for tau function 

Our final case study is focused on the DKP hierarchy. This is a subsystem of 
one of Jimbo and Miwa's integrable hierarchies with with Doc symmetries [55] . 
The master system proposed by Jimbo and Miwa is formulated in terms of tau 
functions t(s, t) with a discrete variable s € Z and a set of continuous variables 
t = {ti,t2, ■ ■ .). These tau functions satisfy an infinite number of Hirota equa- 
tions (or an equivalent set of bilinear equations of the contour integral type), 
which can be divided to three subsets, namely, the equations for T(2s,i)'s, the 
equations for r(2s -I- 1, t) 's and the equations of the coupled type inbetween. The 
Hirota equations of the first and second types are actually identical; it is this 
system of Hirota equations that we call "the DKP hierarchy" . Thus Jimbo and 
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Miwa's master system contains two copies of the DKP hierarchy as subsystems, 
and may be called a two-component DKP hierarchy. 

Since being a relatively less known member of the big family of integrable 
hierarchies, the DKP hierarchy has been rediscovered with a different name such 
as "the coupled KP hierarchy" [51 ISIIIHI US and "the Pfaff lattice" [501 . 
Kac and van de Leur [IS] called Jimbo and Miwa's master system "the charged 
BKP hierarchy" , and formulated a multi-component hierarchy in which the 
BKP, DKP and charged BKP hierarchies and their multi-component analogues 
are unified. 

Let T{s,t) denote the tau function of the DKP hierarchy. Note that it 
amounts to T{2s,t) or r(2s -|- l,t) in Jimbo and Miwa's master system. This 
tau function satisfies the bilinear equation 

^ps-2s -2g5(t-t ,.)^(^/^ ^/ ^ _ 1^ f _ r^-ii) ^ 0. (5.1) 

Here and in the following, the notations are mostly the same as used for the 
KP hierarchy. The integrals on the left hand side are contour integrals along 
a sufficiently large circle |z| = i? (or a formal algebraic operator extracting the 
coefficient of A~^). As in the other cases, this bilinear equation, too, can be 
converted to the Hirota form 

oo 

/i„(-2a)/i„+2s'-2s-i (A)e<"'^*>r(s, t) ■ t(s' - 1, t) 

oo 

+ /in(2a)/i„+2.-2.'-i(-A)e<"^^*>r(s - l,t) • T{s',t) = 0. (5.2) 

Tl = 

For example, if s' = s + 1, terms linear in a„ give the special Hirota equations 

{DiDn - 2hn+l{Dt))T{s,t) ■ T{s,t) 

- 2/i„_3(- A)t(s - 1, t) • r(s + 1, t) = (5.3) 

for n = 3, 4, . . .. 

Note that the left hand side of (j5.3p coincide with that of the special Hirota 
equations (|2.4p of the KP hierarchy. In this respect, the right hand side may 
be thought of as a coupling term among the tau functions r(s, t) with different 
indices; the name "coupled KP hierarchy" originates in such an interpretation. 
This is, however, a rather superficial analogy. Both systems are considerably 
different. For example, as the name "Pfaff lattice" indicates, its relevant equa- 
tions and solutions exhibit a Pfaffian structure [551 [571 155 1 [501 [S^ [M l [S5] 
unlike the determinantal structure of the KP hierarchy. 

Let us mention that the BKP hierarchy is also known to have a Pfaffian 
structure. The difference between the BKP and DKP hierarchies stems from 
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the type of free fermions behind the tau functions [551 SS] • Namely, whereas the 
BKP hierarchy is related to neutral fermions, the DKP hierarchy is formulated 
in terms of charged fermions. The charged fermion system is also used for the 
KP hierarchy; this explains the superficial similarity between the KP and DKP 
hierarchies. Actually, Jimbo and Miwa [55] used three types of free fermion sys- 
tems (charged fermions, neutral fermions and two-component charged fermions) 
to formulate three integrable hierarchies with Doo symmetries. The DKP hier- 
archy is embedded in their first hierarchy formulated by charged fermions. It is 
interesting that their second hierarchy, formulated by neutral fermions, contains 
two copies of the two-component BKP hierarchy. 

5.2 Differential Fay identities 

Two differential Fay identities are obtained from the bilinear equation (jS.ip by 
the following procedure: 

(i) Differentiate the bilinear equation by t'l and specialize s', t' as s' = s + 1 



(ii) Differentiate the bilinear equation by t'l and specialize s' ,t' as s' — s and 



and t' = t+[A-i] + [/^-i]. 



t' = i+[A-i] + [^-i]. 



By (i), the bilinear equation reduces to 




Residue calculus and some algebra yield the equation 



T(,S,f+[A-l] + [/i-^])T(g,f) 
T(,S,t+[A-l])T(s,t +[/.-!]) 



1 T(s + l,t+[A-i] + [/i-i])r(s-l,f) 
AV T(s,t+[A-i])r(s,t + [/i-i]) 




(5.4) 
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In the same way, (ii) leads to the equation 



A - /i t{s, t + [A-i])t(s -l,t~ [/x-i]) 

T(s,t + [A-i]-[M-i])T(s-i,i) 



A-//T(s,t + [A-i])t(s - l,t~ 

r(.,< + [A-i]) 



X + fi- dilog 



r(s-l,t-[M-i])' 



but from the point of view of symmetry, it will be better to shift t as t 
t + This gives the equation 

A2 T(s,t+[/i-i])T(s-l,t+[A-i]) 



X-fi t{s, t + [A-1] + [^l-^)T{s -l,t) 

_ T{s,t+[X~^])Tis-l,t+[fi-^]) 

A-MT(s,t + [A-1] + [^i-^)T{s-^,t) 

(|5.4p was first derived by Adler, Horozov and van Moerbeke in their study 
of the Pfaff lattice . As they pointed out, this equation resembles the 
differential Fay identity (|2.7p of the KP hierarchy; this is another aspect of the 
aforementioned similarity with the KP hierarchy. The only difference is the 
presence of the second term on the left hand side. This is, however, an essential 
difference. 

(|5.5[) has no analogue in the KP hierarchy, etc. A remarkable characteristic 
of this equation and (|5.4[) is that there are two terms, rather than just one, that 
contain the double shift t + [A^^] + [lJ,~^] of t. This apparently small difference 
eventually affects every aspect of this system and its dispersionless analogue. 

Lastly, using the operator D{z), we can rewrite these differential Fay iden- 
tities as follows: 



exp 



((e^(A)_i)(eC(M)_i)iog^) 

- exp((e^(^)+«^^ - l)(e^(^)+^= - l)e-«= logr) 



9i(e^(^)-e^(^))logr , ^, 
T : (5-6) 



A2 

■exp I 



_ _^expf-(e^(^) -e-^Ole''^^^ - l)logr) 



= \ + fi-di (e^(^)+^(^) - e~^^ ) log T. (5.7) 
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5.3 Auxiliary linear equations 

The rather comphcated structure of the differential Fay identities is closely 
related to the fact that an underlying auxiliary linear problem is vector-valued. 
This auxiliary linear problem is substantially the same as the one first discovered 
by Kakei [58] in an inverse scattering formalism, and just a special case of the 
work of kac and van de Leur j45] on multi-component DKP/BKP hierarchies. 
We introduce the wave functions 

T[S,t) 
T{s,t) 



and the duals 



T{s + I,t+[Z-^]) _ 



2s-2 -iit.z) 



-z e 



T{s-I,t+[Z ^]) _^_2s„-at,z) 



T{s,t) 

By suitably shifting s' and s, the bilinear equation ()5.1|) for the tau function 
yields the bilinear equations 

^^^i{s\t\z)n{sAz)+ <f ^^nis',t,z)^,{s,t,z)^0, 

/dz /■ dz 

— *i(s',t',z)^';(s,t,z) + i —■^l{s',t\z)-^2{s.t,z) = 0, 
IiTi J Im 

/dz f dz 

— ■^2{s',t',z)^l{s,t,z)+(j, — *;(s',i',z)*i(s,t,z) = 0, 

for these wave functions. Moreover, these four equations can be cast into the 
matrix form 

rd^f^,{s',t',z) ^l{s',t\z)\ f nis,t,z) n{s,t,z)\_ 

f 21:1 \ *2(s',t',z) ) \ V|/l(s,t,z) ^2{SXZ) ^^-"^ 

This matrix bilinear equation and its building blocks are reminiscent of those 
of the charged two-component KP hierarchy [30j . There is, however, an unneg- 
ligible difference. In the case of the two-component KP hierarchy with charge 
(2s, —2s), the (2, 1) and (1, 2) elements of the matrix wave functions should have 
and z~^*~^ rather than z^''^^ and z^^'^^^ in the foregoing definition of 
the wave functions. Therefore if one interprets this system as a reduction of the 
two-component KP hierarchy, one should remember that this is a very special 
reduction. Actually, Adler, Horozov and van Moerbeke [60] obtained the Pfaff 
lattice in a "deep stratum" of the phase space of the Toda hierarchy. 
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As Kac and van de Leur did in a more general setting [iS], one can derive 
auxiliary linear equations of the form 



(5.9) 



where A„ = An{s,t, di), B,i = Bn{s,t,di),Cn = Cn{s,t,di),Dn = D„(s,t, 5i) 
are differential operator with respect to ti. One can specify the structure of 
these operators in more detail by introducing dressing operators for the wave 
functions. In particular, as Kakei observed [58j . A„ , -B„ , C„ , Z?„ turn out to be 
operators of the form 

An - + o(ar^), Bn - o{d''r^), 
Cn = o(9r'), Dn = -(-ai)" + o(ar^) 

and satisfy the algebraic relations 

An — —Dm Dn = —An, 
^n ~ ~ ^"i 

where A*, B*^,Cn,D*^ denote the formal adjoint of An, -B„, C„, £>„. 



(5.10) 



(5.11) 



5.4 Differential Fay identities and auxiliary linear equa- 
tions 

We now translate the differential Fay identities (|5.4p and (|5.5p to the language 
of the wave functions. As it turns out, each of the differential Fay identities can 
be converted to two different equations (hence altgether four equations) for the 
wave functions. These equations give a generating functional expression of the 
auxiliary linear equations presented above. 

Let us first consider (|5.4p . Shifting t as t — > t — [A^^] — [/i^^], we have the 
equation 

T{s,t-[X-^]-[fi-^])T{s,t) 

T(s,t-[A-i])r(s,*-[/i-i]) 

1 r(. + l,t)T(g-l,f-[A"i]-[/i-i]) 
AV T(s,t-[A-i])r(,s,t-[M-i]) 

= 1 + T di log ' _^ , 

A-^ T{s,t- [n 1]) 

and multiplying this equation by (A— ^)^^''e^*^*'''V(s,t— [/i^^])/T(s, t), we obtain 
the equation 



T(s,t) 

= {di\og^iis,t,\)-di)-9iis,t,fi). (5.12) 
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Similarly, multiplying both hand side of (|5.4D (without shifting t) by (A — 
li)fj.~'^^e~^^*''^'>T{s,t + [fj.~^])/T{s,t) yields the equation 

r{s,t) 

= (5i-9ilogvI/*(s,t,A))vE'^(s,t,Ai). (5.13) 

Note that, unlike the previous cases, there are two terms containing e^^'^'^^. 
They stem from the two terms in (|5.4p containing the double shift t+[A^^] + [/i^^] 
oft. 

We now consider (15. 5|) or, rather its original form 



A^ T(g,t)T(s-l,t+[A-i] - [/i-i]) 
A -/ir(s,t+ [A-i])r(s - l,t- [/i-i]) 

/x^ T(g,t+[A-i]-[M-^])T(s-l,t) 
A - At t(s, t + [A-i])r(s - 1, t - [/^-i]) 

r(s,t + [A-i]) 



A + /X - 9i log 



r(s-l,t-[/i-i])- 



Multiplying this equation by A~^*e"^'^*'^V(s, t + [A"^])/r(s - l,t) yields the 
equation 



^ 'i-lM-e-^(^)*;(., t, A) - ^e-^(^)*^ (. - 1, t, A) 
r(s - l,t) 



- (9i-9ilog^'i(s-l,t,Ai))*J(s-l,t,A). 

Exchanging A and fi and shifting s as s ^ s + 1, we eventually obtain the 
equation 

Ae-^(^)^t (., t, - A-^ + );*^ e-^(^) vl;* (., t, ^) 

r(s,t) 

= (ailog*i(s,t,A)-ai)vI/t(s,t,^), (5.14) 

which has the same structure as (|5.12p . Similarly, multiplying the same equation 
as above by /z^*'~^e^(*'''V(s — l,t — [^~^])/T(s,t), we obtain the equation 

Xe^('H,{s, t, m) - A-i ^i^^^e^(^) *i(., t, ^,) 

T[S,t} 

= (91 -9ilog«'^(s,t,A))^'2(s,t,^), (5.15) 

which can be compared with (|5.13p . 

These results (j5.12p - ()5.15p show that the vector-valued wave functions 

$1 \ / *i(s,t,Ai) \ / **(,,<, ^) 

$2 / I *2(s,t,Ai) J ' I n{s,t,fi) 
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satisfy the same linear equations 

T{s,t) 

Ae^(^)$2 - A-i^^^f^e^(^)<i>i = (^i - logM/*(s,*, A))$2. 

T{s,t) 

Expanded in powers of A, these equations give rise to an infinite number of 
Unear equations of the form 



T{s,t} 
T(S,t) 

for n = 2, 3, . . ., where Vn and u„ are defined by 

Vn = dihn-2{-dt)logT{s,t), Un = 5i /l„_2 (9t ) log t(s, t) . 

The lowest equations (for n — 2) read 
2 T{s,t) 

lidf + d2)<S>2 - ^^^7^<fi = -(a?logr(s,t))cI>2, 
2 T{s,t) 

which one can rewrite as 



(5.16) 



2li^ af-25?logr(s,<) j V * 



^2 



This is exactly the lowest equation of the auxihary hnear problem (|5.9p . Higher 
equations can be derived recursively by the same procedure as in the case of the 
KP hierarchy etc. 

5.5 Quasi-classical limit — encounter with difficulty 

Quasi-classical limit can be achieved by assuming the ansatz 

Tn{s,t) = exp(h~^F{s,t) + 0{h~^)'j (5.17) 

for the rescaled tau function Tfi{s,t) = T{h,h^^s,h~^t). The differential Fay 
identities, written in the form of l|5.6p and (|5.7p . yield the dispersionless Hirota 
equations 

gi3(A)r>(M)F _ y^-2^-2^{D(\)+^s)(D{^.)+^,)F ^ ^ _ di{D{X) ~ D{^i))F ^ ^^^^^^ 

X — fi ' 
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-^^ ,-_D(A)(D(p)+a,)F _ _J^-{D(\)+d^,)D{p)F 



A — /i A — fj, 

= \ + fi-di{D{X)+D{fi)+ds)F (5.19) 

for the F function. The first equation obviously resembles its counterpart p.23p 
for the KP hierarchy; the only difference is the presence of the second term on 
the left hand side. On the other hand, the second equation seems to have no 
analogue in the previous cases. 

Let us examine these equations in more detail. It is convenient (and natural) 
to introduce the S function 

S{z) = £,{t, z) + 2s log z - D{z)F 

and rewrite the two equations as 

^D{\)D(p)F(^^ _ ga>-a,s(A)-a,s(M)) ^ diSjX) - diSjfi) ^ ^g) 

A — fi 

-D(X)D{f^)F 

^ (^e^^sw _ g9,s(M)) ^ Q^s{X) + diSifi) - dsdiF. (5.21) 

A — 

One can eliminate e^'^^^^'^^^^^ by multiplying both hand sides of these equations. 
This yields the relation 

which one can rewrite as 

(5i5(A))2 - {d,d,F)idiS{X)) - e^=^(^) - e^=^-^=^(^) 

Since A and /i are separated to the left hand side and the right hand side, 
both hand sides of this equation are actually independent of A and /i. One can 
determine their value by letting A, — > oo. We thus obtain the equation 

i^lS{^l)f - ididsF){diS{n)) - e^=^(^) - e^'^-^=^(^) 

= ~2d^,F + ^d2dsF - i(ai9,F)2 (5.22) 

as a consequence of (|5.20l) and (|5.2ip . Actually, this procedure can be reversed. 
Namely, assuming that (|5.20p holds, one can easily recover (|5.2ip from (|5.22p . 
Thus under ((On|) . ((OT|) and ((O^ are equivalent. 

(|5.22p may be thought of as a quadratic constraint to (|5.20p . This constraint 
can be solved for either 9i5(/i) or e^''^'^-', which becomes a non-polynomial 
function of the other. For example, solving it for e^='^*^''' gives 

^OsSi^.) ^ 1 |^p(^)2 _ (^dAF)p{p) -I + Vd), (5.23) 
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where 



pifi) = d^Sifi), D = {p{i,f - idid,F)p{ti) - /) - 4e^>, 

and / denotes the right hand side of (|5.22p . This enables one to ehminate e^='^(^) 
and e^'"^('') from dOOl)- The outcome is an equation for L>(A)L>(^)F, diS{X) 
and diS{^). One can rewrite it in a form that amounts to (|2.3ip . Thus S'(/i) 
anyhow turns out to satisfy a set of Hamilton- Jacobi equations of the form 

dnSifi) = Hn{diS{fi)). (5.24) 

A serious problem shows up here. Namely, unlike the cases of the KP, Toda 
and BKP hierarchies, the right hand side of these Hamilton-Jacobi equations 
is not a polynomial function of diS{fi) but a fairly complicated (though alge- 
braic) irrational function thereof. Roughly speaking, Hn{p) is a combination of 
polynomials and the square root of a quartic polynomial, namely, 

Hn{p) = Pn{p) + Qn{p)^, (5.25) 

where P„ (p) and Qn (p) are polynomials in p, and D is given by 
D^ip'~{d,d,F)p-f)^-4e^'^^. 

This new phenomena, hinting at the relevance of an elliptic curve, is closely 
related to the fact that the auxiliary linear problem (|5.9[) is essentially vector- 
valued. Because of this complicated structure of the Hamiltonians Hn{p), the 
question of integrability of the underlying dispersionless system has not been 
resolved. 

These results could have been derived by the quasi-classical limit of the aux- 
iliary linear equations. Such an approach, however, is technically more subtle, 
because the auxiliary linear problem has matrix coefhcients. The foregoing ap- 
proach based on the differential Fay identities is obviously simpler and more 
reliable. This is a demonstration of the advantage of using differential Fay 
identities. 

6 Concluding remarks 

One of the most remarkable results of our case studies is that differential (or 
difference) Fay identities are the auxiliary linear problem in disguise. This 
interpretation clarifies the meaning of the mysterious linear equations (|2.19p 
of Sato and Sato ^] and its generalization to other integrable hierarchies. 
The status of differential Fay identities is thus parallel to that of dispersionless 
Hirota equations. We can use them as a new foundation of dispersive integrable 
hierarchies. In this new framework, we might be able to go beyond dispersionless 
(or quasi-classical) limit to study higher orders of fi-expansion. 

Though the case studies in this paper are limited to the KP, Toda, BKP and 
DKP hierarchies, these results can be (and have been) generalized to some other 
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cases. A particularly successful case [55] is the charged multi-component KP 
hierarchy and its dispersionless analogue (which is identified to be the genus- 
zero universal Whitham hierarchy |10|V In a sense, this is a multidimensional 
generalization of the Toda hierarchy as well. The results on the KP and Toda 
hierarchies can be fully generalized to this case. 

Among the cases already studied, the DKP hierarchy and its dispersionless 
analogue are still posing many open problems. The technical complexity in 
this case stems from the fact that the auxiliary linear problem is essentially 
vector- valued. In contrast, the charged multi-component KP hierarchy has ac- 
tually a scalar-valued auxiliary linear problem, by which one can formulate its 
dispersionless analogue as the genus-zero universal Whitham hierarchy. Our pre- 
liminary consideration in this paper indicates that the dispersionless analogue 
of the DKP hierarchy might be related to the universal Whitham hierarchy of 
genus one [1^ or fall into one of examples of "quasi-classical deformations of 
algebraic curves" proposed by Kodama, Konopelchenko and Marninez Alonso 
[57] . The Toda versions of the DKP hierarchy [55] [53] will inherit a similar 
problem. These issues deserve to be further studied. 

To conclude this paper, let us mention a geometric interpretation of disper- 
sionless Hirota equations. This interpretation is inspired by the recent work 
of Krichever, Marshakov and Zabrodin [70] (see also the work of Carroll and 
Kodama [12] ) . As the recent work of Eynard and Orantin [7T] on random ma- 
trices indicates, such a geometric point of view will be indispensable when one 
pursues higher orders of fi-expansion of log . For simplicity, we now consider 
the case of the dispersionless KP hierarchy, but it is rather straightforward to 
generalize this interpretation to the dispersionless Toda hierarchy and the uni- 
versal Whitham hierarchy of genus zero. A clue is to take the limit as ^ — > A in 
the dispersionless Hirota equation 



e 



d(a)d(m)f ^ p(A) -p(/i) 
A — /I 



This yields the equation 

by which one can rewrite the dispersionless Hirota equation as 
Remarkably, the quantity 



E{\,^i) 



showing up on the right hand side is nothing but a local expression of the prime 
form [32l [72] on the Riemann sphere in a neighborhood of z ^ oo with local 



44 



coordinate p{z) ^. We can thereby rewrite (|6.ip as 

Ei\, = (A - exp(-ip(A) - . (6.2) 

This coincides with a formula presented by Krichever et al. [70] for the case 
of nonzero genera. Moreover, taking the logarithmic second derivative of both 
hand sides with respect to A and fi, we obtain an expression of the Bergmann 
kernel B{X, /i) = dxd^ \ogE{X, /x) [321E2]: 

B{\, A*) = + D'{X)D'{fi)F. (6.3) 

Remarkably, this formula resembles a well known result on the two-point loop 
correlation function (which corresponds to D' {X)D' {fi)F in this formula) of the 
Hermitian random matrix model in the large- limit. This seems to indicates, 
as Eynard and Orantin [71' observed in the case of random matrices, that the 
Bergmann kernel is one of fundamental building blocks of ?i-expansion of logrn- 
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